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Abstract 

The odd A proton drip line nucleus ^^F has been studied in a semi-microscopic 
model which couples proton quasiparticle motion to the vibrational motion of 
the neighbouring even-even ^^Ne core. The experimentally observed low lying 
excitation spectrum, electric quadrupole moment, magnetic dipole moment, 
B{E1) and B{E2) values have been fairly well reproducd. The calculated rms 
radius of the first excited state is well reproduced and is found to be larger 
than that of the ground state which agrees with the experimental observation. 

PACS number(s): 21.60 Ev, 27.20 +j 
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I. INTRODUCTION 



With the advent of radioactive ion beam facihties, new experimental data concerning 
exotic nuclei are becoming available. Some of the exotic nuclei show halo structures which 
has generated a lot of interest, in recent times, among various theoretical and experimental 
groups. Experimental investigations have shown the enhancement of sub-barrier fusion 
cross-section in reactions involving exotic nuclei which is thought to be related to the 
sizes and structure of such nuclei. Thus the understanding of the structure of such nuclei is 
of utmost importance to the study of various phenomena involving these nuclei. 

It is well known that the very low binding energy of the last unpaired nucleon is respon- 
sible for the large matter radii of one-nucleon halo nuclei. A number of one-neutron halo 
nuclei have been identified and studied in recent times. In a similar way, one-proton halo 
nuclei ought to exist except for the fact that the proton wave function, unlike the neutron 
wave function, is squeezed due to the Coulomb barrier. The nucleus ^^F may be a good 
candidate for a proton halo nucleus since the binding energy of the last proton is only 0.60 
MeV. Experimental studies 0-|5| have revealed that the spin-parity of the ground state of 
^^F is I and has spectroscopic strength equal to 0.93. The first excited state is known to be 
i^, with a dominant 2si/2 single particle (s.p.) configuration having binding energy equal to 
0.1 MeV. However, the root mean square (rms) radius of the ground state is 3.78 fm whereas 
that of the first excited state, at 500 keV excitation, is observed to be 5.33 fm Consid- 
ering the ground state and the first excited state to be pure single particle configurations 
of ld5/2 and 2si/2, respectively, the rms radii of the last proton, using Harmonic oscillator 
wave functions, turn out to be the same, which is contrary to the experimental findings. 
This indicates that coupling of the odd proton motion and that of the core may modify the 
wave function of ^''F sufficiently to affect the rms radius of the ground state and that of the 
first excited state. The neighbouring even-even cores for this nucleus are ^^O and ^^Ne. The 
first, ^^O, is a doubly closed shell nucleus without any vibrational structure, whereas the 
second, ^^Ne, shows a well developed vibrational structure with the first 2"*" state at 1.88 
MeV [0. As such, the behaviour of the core changes abruptly if one goes from O to Ne. 
Therefore, the nature of the ^''F nucleus is in between the two extremes and the structure 
of this nucleus needs to be explored. Furnstahl and Price 0, through a relativistic Hartree 
calculation, found the rms charge radius to be 2.76 fm and the magnetic dipole moment 
to be 4.87 fiN- On the other hand, the density functional method of Lombard provided 
the total binding energy (B.E. = 126.6 MeV). However, the detailed structure information 
regarding the excitation mechanism is lacking at present. 

We have shown earlier that the structures of the one-neutron halo nuclei ^^Be and 
^^O could be explained satisfactorily through the coupling of collective vibrations of the 
respective even-even cores with the motion of the odd valence neutron. In the present 
work, we have shown that similar idea can be applied in explaining the important features 
of the excitation spectrum of the proton rich exotic nucleus ^^F by coupling single proton 
qusiparticle motion to the collective vibrational motion of the even-even ^^Ne core. In this 
paper, we present the results of our calculations for the studies of the low lying excitation 
spectrum and the electromagnetic (EM) properties of the ^^F nucleus. The model used in 
the present calculations is discussed briefiy in section II. The method of calculation is given 
in section III. Section IV contains the results and discussion and finally the conclusions are 
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given in section V. 



II. MODEL 



The model used in the present calculations is given in our earlier work [|10],|rT| and here 
we are giving it just for the sake of completeness. The total Hamiltonian for a coupled 
system of a quasiparticle and collective excitation of the core may be written as 

H = Hcore + Hs,p, + Hint (l) 

where Hcore describes the harmonic collective vibrations of the even-even core. Hg.p, describes 
the motion of the valence nucleon (-hole) in an effective potential. The basis states, which 
are eigen functions of Hq = Hcore + Hs,p,, can be written as {[N2V2L2, NsVsL3]^{nll/2y}-^^ , 
where {NxV\Lx) is the totally symmetric state of N\ phonons of multipolarity A coupled 
to angular momentum Lx and Vx represents all the additional quantum numbers necessary 
for complete specification of the state; L2 and L3 are coupled to the angular momentum L 
of the core which is then coupled to j of the quasiparticle to form the resultant J. (n/1/2) 
represents a quasiparticle state having radial quantum number n, orbital qauntum number 
/ and spin 1/2. The eigen values of Hg are the sum of the single particle energies and the 
core energies. The core particle interaction is given by 

23 X 

H.nt = E Kxir) E cyl,Yx,iep, 0,), (2) 

A fj,=—X 

Yx^ are the spherical harmonics corresponding to the particle coordinates. The collective 
coordinates ax/i can be expressed as a linear combination of the phonon creation and de- 
struction operators b respectively, for the core. The matrix elements of the interaction 
Hamiltonian, after taking care of the pairing effects, are given by 

< {[N2V2L2,N3V3L,]'^{nll/2yy\H,^t\ 

{[Kv^L',, N;,viL',fin'n/2y'y'' > 

= E^A < n/||[(^)'\]^||n7' > (-1)^-V2+A 

A 

x^x{[L][L'][j][/]}^x 

L+L'+La+L' 





A-l 



1)L2+L3^ NxLx\\bt\\NyLx^ > 



+Sn',n,^, X (-l)^'.+^3 < N',L',\\bt\\NxLx > 

In the above, X2 and X3 are the quadrupole and octupole interaction strengths respectively, 
and 7] = 2 when A = 3 and vice versa. Xa's are given by 
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X^ = K^J2{\ + 1\—^, (3) 



where Kx are the average values of the radial integrals. The quantities V^ 's and t/^'s are the 
occupation and non-occupation probabilities of the respective single particle states. Diago- 
nalising the Hamiltonian of the coupled system in the basis space of Hq., one obtains the eigen 
values and the eigen vectors. The eigen vectors, thus obtained, are used to calculate the 
static electric quadrupole and magnetic dipole moments, the electric and the magnetic mul- 
tipole transitions. The expressions for the static magnetic dipole and the electric multipole 



operators for the coupled system are |12 



0(M1, ^l) = (3/47r)^/' {gil^ + g^s, + gnR) , (4) 

0(i?A,/i) = |-Zei?o^ (^)^ {hi, + (-)^&A,-,) + (e, + ^) Y,,,{e,<j^) (5) 

where Z, A and Rq are the charge, mass and the radius of the core, respectively. The quantity 
Cx refers to the stiffness constant of the vibrating core. The magnetic dipole moment /i and 
the electric quadrupole moment (QM) Q of a state | J, M, a > are defined as 

Air 

fi=< J,M = J,a\ {—) O{Ml,0) \J,M = J,a> (6) 
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and 



IQtt 1/2 

Q =< J,M = J,a\ {—) O{E2,0) \J,M = J,a> . (7) 



The reduced eletromagnetic transitions are given as 

B{EX : J, ^ Jf) =<Jf\ \ 0{EX) II >2 {2Ji + 1)"^ (8) 

The detailed expressions, after carrying out angular momentum algebra, are given in 
Heyde and Brussaard [1^ except for a multiplicative factor in the particle part invoving 



Vj's and U/s. The rms radius of a state is obtained by taking the expectation value of the 
operator r^, in the respective state. The core contribution to the rms radius is either taken 

from the experiment, if available. Otherwise, it is taken as Rq^ = {1.2A^^^)'^ . 



III. METHOD OF CALCULATION 



There are several parameters in this calculation viz. (i) the quasiparticle energies ej, (ii) 
the Uj{Vj) factors and (iii) the interaction coupling strengths X2 and X3. The vibrational 
energies for the i^Ne core are taken from the experimental data 0. The experimental 
excitation spectrum of i^Ne shows the behavior of a vibrator. The first 2"*" excited state 
occurring at 1.887 MeV can be identified as one quadrupole phonon excited state. The 
second 0"*" , 2+ states and the first 4"'" state occurring at 3.576 MeV, 3.616 MeV and 3.376 
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MeV, respectively, can be taken as the members of the multiplet of two quadrupole phonon 
excitations. The calculations have been restricted to only two quadrupole phonon space as 
the calculated spectra have been found to be insensitive to the inclusion of octupole phonons. 
Therefore, X3 is kept equal to zero. The single particle space, in the present calculations, 
is confined to 1^5/2, 2si/2 and 1^3/2 orbitals for the positive parity states and lp3/2, lpi/2 
orbitals for the negative parity states. The Uj{Vj) factors for the 1^5/2, 2si/2 and 1^3/2 single 
particle states are available from the {a, t) and {^He, d) reaction studies |T^,^ . However, 
the experimental spectroscopic information for the IP3/2 and lpi/2 single particle states is 
not very clear. From the (rf,^ He) reaction, Firestone et al. |]T3| overestimated the strengths 
for the lp3/2 and lpi/2 single particle states by 30 — 40% and, to the best of our knowledge, 
no other data are available. Therefore, we have taken them as free parameters. The UjiVj) 
factors used in the present calculation to reproduce the experimental data are given in Table 
I. The relative quasiparticle energies, e^, with respect to the lowest single particle state in 
the conventional shell model basis and the interaction coupling strength X2 are treated as 
free parameters. These are varied to fit the experimental level scheme and the spectroscopic 
factors. The set of best fit parameters used in the present calculations are given in Table |. 

The wave function of a state with angular momentum J with z-component M and energy 
i?" (a distinguishes between states of same spin and parity) is given by 

|E",JM>= Co^{N2L2,j,J)\N2L2,j;JM>, (9) 

and the spectroscopic factor is defined as 

5"(/,j = J) = t/2|C„(00,j,J)p. (10) 



IV. RESULTS AND DISCUSSION 



A. Even parity states 



For the even parity states, the odd quasiparticle motion was restricted to be in the Ids, 

2s 1 and Ida orbitals. The collective states of ^^Ne, used in the present calculation, to which 

2 2 

the quasiparticle couples, were taken from the experimental data 0. The calculated low 
lying levels and the spectroscopic factors are shown in Fig.l, along with the experimental 
data, [0 , for the sake of comparison. It is seen from this figure that the present calculation 



is able to reproduce the positive parity level spectrum quite well. The present calculation 
predicts two weak and states, not observed experimentally, at ~2 MeV excitation 
with the spectroscopic strengths 0.05 and 0.06, respectively (not shown in Fig. 1). . The 

Q _|_ 

calculated spectrum is found to have three | states lying between 4.50 and 4.75 MeV 
excitation with single particle strengths equal to 0.05, 0.51 and 0.07, respectively. The state 
at 4.63 MeV with a spectroscopic strength of 0.51 may correspond to the experimentally 
observed |^ state at 5 MeV excitation with spectroscopic strength 0.54. Whereas, the 
calculated state at 4.75 MeV may correspond to the observed state at 5.82 MeV. The 
calculated I"*" and |^ states at 6.83 MeV and 6.84 MeV may correspond to the observed 
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levels at 6.56 MeV and 7.36 MeV, respectively. It is conjectured that the experimentally 
observed I"*" state at 8.07 MeV may be associated with the calculated state at 8.20 MeV. 
There is an overall good agreement between the calculated and the experimental level 

Q _|_ 

energies and the spectroscopic factors. However, almost all the calculated | states are 
depressed in energy by 500 keV as compared to the correponding experimentally observed 
states. The observed Ids single particle strength is found to be almost totally concentrated 

on the Ij*" state though the present calculation predicts a |^ state at 2 MeV with a very 
small spectroscopic strength. A similar situation is observed with the Ida single particle 

2 

strength distributions. The available Irfs single particle strength is found to be fragmented 
amongst three states lying between 4.50 MeV and 4.75 MeV, the major part of strength 
being concentrated on the state at 4.63 MeV excitation. 



B. Odd parity states 

For generation of the odd parity states, the single particle space considered in the present 
calculation is spanned by Ipi and Ipa states. The quadrupole interaction coupling strength 
X2 was kept identical to that used in the case of the even parity states and the single particle 
energies eip^ and eip^ were treated as free parameters. The calculated level spectrum of the 

low lying odd parity states is shown in Fig. 2 alongwith the experimental one [Q. It 
is seen from this figure that there is an overall good agreement between the calculated 
and the experimental spectra except for the | states which are overpredicted by ~ 1 

MeV. The calculation predicts that almost the total available Ipi single particle strength 

2 

is concentrated on the | state at 3.186 MeV. The Ipa single particle strength is predicted 

^ 2 

to be fragmented into three | states lying at 4.74 MeV, 5.20 MeV and 7.33 MeV with 
major part of the strength going to the state at 4.74 MeV. The structure of the calculated 
I state at 4.86 MeV is found to be mainly consisting of a 2+ collective state coupled to the 
Ipi single particle orbital. Similar are the structures of the | and | states. However, no 
experimental data are available for the spectroscopic strengths for the odd parity states of 
^^F to compare with. 



C. RMS radii, EM moments and Transition probabilities 



The stringent test of the wave functions obtained by diagonalising the Hamiltonian is to 
see whether they are able to reproduce the electromagnetic moments and the multipole elec- 
tromagnetic transition probabilities. However, in the case of the exotic nuclei an important 
aspect of structure study is the reproduction of the experimentally observed rms radius of 
proton or neutron distribution which provides a measure of the proton or neutron halo. The 
rms radii of the ground state and the excited states of ^''F contains two parts - (i) from the 
core and (ii) from the extra quasiparticle. In order to evaluate the proton rms radii of these 
two states, initially, we performed our calculation with the Harmonic Oscillator (H.O.) wave 
functions to extract the contribution from the quasiparticle part. The calculated values of 
the proton rms radius for the |^ ground state and the first excited state were found to be 



6 



4.01 fm and 4.15 fm, respectively. Though the experimental trend for the rms radii for these 
two states is reproduced, the calculated values are not in agreement with the experimentally 
observed values of 3.78 fm and 5.33 fm, respectively 0. This discrepancy may be attributed 
to the fact that the H.O. wave functions do not carry any effect of the low binding enegy 
of the last valence nucleon. This behaviour of the H.O. wave functions, through the radial 
integrals, consequently affects the calculated values of the electric quadrupole moment and 
the electromagnetic transition rates (Eq.^. 

To overcome this difficulty, we have used a generalised average one-body potential of 
Woods-Saxon type |]T^ to generate the s.p. wave functions. Keeping all the parameters 
fixed, the depth Vq and the radius parameter ro were varied to reproduce the experimentally 
observed binding energies, 0.600 MeV and 0.105 MeV of the Ids and 2si s.p. states, 
respectively. The parameters of the single particle potential are given in Table y. These 
W-S wave functions could be used for calculating the rms radius and the electromagnetic 
properties of the states of ^^F but for the wave function of the Ida state, which is unbound. 

In order to circumvent this problem to get the wave function of the unbound Ids s.p. 

2 

state, we developed the W-S wave function in the H.O. basis (WSHO). This was done by 
digonalising the W-S potential (Table 0) for each (/, j) state in H.O. basis states {n,l,j) 
with nodal quantum number n=10 and properly picking up the Ids, 2s i, Ida, Ips and Ipi 
states. The < > values for the even parity s.p. states, thus developed, are given in Table 
rri. From the table it is observed that < > values calculated with the W-S wave function 



in the H.O. basis are smaller than those calculated with the W-S wave functions and quite 
larger than those calculated with the H.O. wave functions. The calculated values of the rms 
radii of ^ ''F using the quasiparticle wave functions WSHO are given in Table |IV|. For the first 
excited state, the calculated value of the rms radius agrees quite well with the experimental 
one. However, for the ground state the calculated value of the rms radius overestimates the 
corresponding experimental value. 

The same quasiparticle wave functions, as used in the calculation of the rms radii, have 
been used to calculate the electric quadrupole moment of the ground state and the EM 
transition probabilities among the various states of ^^F obtained earlier. Effective proton 
charge, (ep)gjj = 1-5 e and effective charge of the core, Z^fj = 0.2 e have been used in the 
present calculation (e being the charge of proton). However, we made a rough estimate of 
the value of Z^ff (defined as Z{hu)2l1C2)^^'^) from the value of the quadrupole interaction 
strength parameter X2 ( see Eq. |^ ) used earlier. Taking value of Uuj^-, the quadrupole phonon 
energy ^ 2 MeV, and of < -ft" >, the average radial integral ~ 40 MeV, the calculated value 
of Ze// is found to be equal to 0.20 e, almost same as the one used in the calculations of the 
QM and the EM transition probabilities. The calculated values of the QM of the ground 
state and the reduced transition probabilities B{E2) and B{E1) are given in Tables |^ and 
|V| along with their corresponding experimentally measured values |T6|. For evaluation of the 
magnetic dipole moment, the values of gs, gi and gn were 5.585, 1.0 and 1.0, respectively 
The calculated values of the electric quadrupole moment and the magnetic dipole 



moment are 0.108 eb and 4.708 /xat, respectively and they compare very well with the 
experimental values of 0.10 ±0.02 eb and 4.72130 ± 0.00025 fiN, respectively (see Table 
fV]). The calculation is able to reproduce well the observed reduced transition probabilities, 
B[E2) and B[E1) (Table 0). However, the calculated B[E1) value between the |^ and |^ 
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states underestimate the experimental value by an order of magnitude. 

V. SUMMARY AND CONCLUSION 

The low lying excitation spectrum of the proton drip line nucleus ^''F has been studied in 
the quasiparticle-core coupling model in which the motion of the odd quasiparticle is coupled 
with the neighbouring even-even vibrating core. The core excitation energies, which take 
into account all the correlations arising out of the many particle core system, were taken from 
the experiment. The relative single particle energies and the quadrupole interaction strength 
were treated as parameters. The calculated energies and the spectroscopic factors for both 
the even and the odd parity states agree quite well with the corresponding experimental 
results. The calculated values of the rms radii of the ground state and the first excited 
state using single particle H.O. wave functions follow the experimental trend but are grossly 
underestimated (particularly the rms radius of the first excited state). However, the value of 
the rms radius of the first excited state is well reproduced by using the WSHO single particle 
wave functions. The calculated values of the electric quadrupole moment of the ground state 
and the reduced transition probabilities using the WSHO wave functions compare quite 
well with the corresponding experimental values. The calculated value of the ground state 
magnetic dipole moment agrees very well with that of the experimental one. 

In conclusion, it may be inferred that within a simple quasiparticle-vibration coupling 
model, one can very well reproduce the experimental observables associated with the struc- 
ture of the light proton drip line nucleus ^''F without going into sophisticated calculations. 
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FIG. 1. Calculated and experimental even parity levels of ^^F. The excitation energies, 
spin-parities and spectroscopic factors of the levels are shown. 
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FIG. 2. Same as in fig. 1 for odd parity levels of ^^F. 
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TABLES 



TABLE L Parameters of the Calculation 


Ids 


2si 


iT,- (MeV) X2 
Ids Ips Ipi Ids 2s 1 Ids Ips Ipi 

9 9 999 9 99 


0.985 


0.866 


0.800 0.142 0.313 0.0 0.550 4.500 4.800 3.00 0.5 


TABLE IL Parameters of Woods-Saxon s.p. potential 


^0 




ro ao Vs rg 


49.091 




1.404 0.685 23.937 0.819 1.060 


TABLE in. Mean square radius for single particle orbitals 


< > 


W-S 


wave fuction in H.O. Basis W.S wave function*^ H.O. wave function 


Ids 

2 
2 

Ids 

2 




14.154 16.012 8.819 
19.342 29.925 8.819 
17.690 8.819 


'^Parameters of Table 1 






TABLE IV. Electromagnetic moments and rms radius 


J"" 


Quadrupole Moment (eb) Magnetic Moment (/ijv) rms radius (fm) 
Calc. Expt. Calc. Expt. Calc. Expt. 


5 + 
2 

1 + 
2 


0.108 


0.10±0.02 4.7080 4.72130±0.00025 4.83 3.78 
2.7706 5.34 5.33 


TABLE V. Reduced Electric transitions probabilities B{E2) and B{E1) 




TTT 

-4 


^i(MeV) Ef{Me\) Mult. Calc. (W.u.) Expt. (W.u.) 


1 + 

2 

1- 
2 

5- 
2 


5 + 
2 

1 + 

2 

5 + 
2 


0.517 0.0 E2 21.9 25.0±0.5 
3.087 0.517 El 0.2x10"^ (1.5±0.3)xl0-3 
4.848 0.0 El 6.1x10-^ (4.3±0.8)xl0-3 
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